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ABSTRACT 


During  late  1964,  a  joint  Canadian  and  United  States  Air  Force  ex¬ 
periment  was  carried  out  at  Cold  Lake,  Alberta,  Canada,  to  examine  multi- 
path  propagation  effects  on  various  types  of  radio  signals  A  relatively 
simple,  vertical  reflection  interferometer  was  used  to  record  the  received 
signal  variations  from  an  airborne  radio  source.  Direct  meteorological 
measurements  were  carried  out  along  the  paths  of  the  radio  signals  within 
the  same  time  frame  This  paper  is  concerned  with  the  radiometric  results 
which  relate  the  radio  signal  fading  pattern,  due  to  multipath  to  effectivt 
earth  radii  values  Using  the  meteorological  measurements  independent 
ray  tracing  calculations  demonstrate  that  the  interferometer  provides 
meaningful  measurements  of  the  effective  ray  path  bending  under  various 
propagation  conditions.  Based  on  these  observations  the  following  paper 
develops  a  method  to  generate  the  effective  radio  refractivity  profile  from 
direct,  interferometer  radiometric  measurements  An  analysis  of  the 
accuracy  requirements  indicates  that  the  method  is  within  the  state-of-the- 
art  and  can  be  modified  in  a  simple  way,  using  two  frequencies,  to  calculate 
the  effective  height  of  the  receiver  relative  to  the  reflecting  surface.  This 
latter  ability  is  particularly  important  for  measurements  over  the  ocean. 
The  technique  appears  to  provide  a  direct  and  simple  method  to  measure 
effective  radio  propagation  conditions  in  real-time.  For  many  applications 
this  approach  has  significant  advantages  since  it  reduces  the  usual  require¬ 
ment  to  take  extensive  airborne  meteorological  measurements 
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1. 


INTRODUCTION 


Radio  ray  propagation  conditions  are  usually  determined  by  direct 
meteorological  measurements  through  the  use  of  balloon -borne  instrumen¬ 
tation  and/or  instrumented  aircraft.  This  approach  has  two  disadvantages 

In  the  first  case,  if  the  propagation  conditions  are  to  be  defined  in 
detail  the  measurements  can  require  the  use  of  several  airborne  facilities  to 
counteract  the  effect  of  spatial  and  temporal  variations  This  approach  leads 
to  a  large  amount  of  data  with  the  usual  difficulty  in  processing  them  in  a  us eful 
time  period 

Secondly,  the  effect  of  small  variations  of  the  radio  refractivity 
structure  on  radio  propagation  parameters,  such  as  range  errors  cannot  be 
determined  without  additional  ray  tracing  analysis 

The  following  presentation  suggests  a  method  to  determine  the  ef¬ 
fective  radio  refractivity  model  by  direct  radio  propagation  measurements 
The  concept  was  generated  by  results  obtained  from  a  comprehensive  field 
experiment  wherein  a  vertical  reflection  interferometer  was  used  to  determine 
the  angle  of  arrival  of  radio  signals  In  the  final  analysis  a  comparison  was 
made  between  the  "effective  earth  radii’1  calculated  from  the  radio  inter¬ 
ferometer  data  and  by  a  separate  calculation  using  direct  meteorological 
measurements  together  with  ray  tracing  analyses 

The  agreement  between  the  comparative  data  was  significantly  good 
and  the  idea  was  then  germinated  as  to  the  possible  direct  use  of  radio  inter¬ 
ferometer  data  to  measure  the  effective  propagation  conditions  without  direct 
meteorological  support  The  final  projection  was  to  determine  if  an  "effective" 
radio  refractivity  profile  could  be  generated  using  the  interferometer  earth 
radii  data  The  following  analyses  will  demonstrate  that  this  is  analytically 
and  technically  feasible 
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THE  DETERMINATION  OF  EFFECTIVE  EARTH  RADII  FROM  RADIO 
INTERFEROMETER  MEASUREMENTS 


In  the  reflection  interferometer  technique  the  signal  from  the  radio 
source  arrives  at  the  receiving  antenna  via  two  paths  (Figure  1),  i.e.,  by 
the  "direct  ray"  (Si)  and  by  the  "reflected  ray"  (Ss).  S1  and  Ss  combine 
vectorially  according  to  their  phase  differences  and  also  their  relative 
amplitudes 

The  phase  difference  between  the  two  signals  is  a  function  of  the  path 
length  difference  of  the  two  rays  in  terms  of  wavelengths  (X)  of  the  trans¬ 
mission  frequency  and  also  of  the  phase  change  in  the  reflected  ray  at  the 
surface.  For  water  and  small  angles  of  incidence  the  phase  change  on  re¬ 
flection  is  essentially  180  degrees  Therefore,  the  two  signals  will  cancel 
at  the  receiver  when  the  electrical  path  length  differences  are  integral 
multiples  of  the  transmission  wavelength  Then  the  received  power  will 
oscillate  between  a  maximum  and  min  mum  level  given  by 

(1  -  p)  P  <  Pr  <  (1  +  P)  P  (1) 

where  P  is  the  power  received  along  the  direct  path,  S1?  and  pis  the  magni¬ 
tude  of  the  reflection  coefficient.  As  the  source  moves  toward  the  receiver 
the  signal  will  be  a  minimum  when  the  path  length  difference,  AR,  is 

AR  =  nX  (2) 

wher e  n  -  1 ,  2,  3,  etc 


Figure  2  shows  a  typical  received  signal  power  pattern  as  the  radio 
source  varies  in  range  with  respect  to  the  receiver  and  maintains  a  constant 
height  from  the  reflecting  surface 


Referring  to  Figure  3,  the  electrical  ranges  ID,  measured  over  the 
direct  path,  which  produce  minima  (fades)  is  given  by 


ID  n  - 


2b!1  hp1 

nX 


(3) 


The  true  heights  of  the  receiving  antenna,  hls  and  source,  hg,  relative 
to  a  curved  earth  are  AA1  and  BB';  respectively  As  shown  in  Appendix  I, 
the  range  IDn  (Equation  (3))  can  be  determined  in  terms  of  ha,  h2,  and  the 
earth  radius,  A 
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FIGURE  I  VERTICAL  INTERFEROMETER  STRAIGHT  RAY  GEOMETRY 
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FIGURE  2  VERTICAL  INTERFEROMETER  INTERFERENCE  PATTERN 
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FIGURE  3  GEOMETRY  OF  FLAT  EARTH  APPROXIMATION 
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In  practice  the  rays  do  not  travel  in  straight  lines,  however,  it  was 
shown  by  Schelleng,  Burrows  and  FerrelA  that  by  increasing  the  earth  radius 
to  some  value  Aes  the  curved  path  geometry  can  be  represented  by  straight 
rays,  preserving  the  initial  elevation  angle  of  the  ray,  0O,  the  electrical  path 
range,  R,  and  the  height  of  the  sourc  hg?  relative  to  this  new  earth  surface 

The  new  "effective  earth  radius",  Ae,  is  calculated  from 


Ae  ~ 


A  n0 


nr 


10 


g-  dN /dh  cos  0 


(4) 


where 


A 

no 

dN/dh 

9 


is  the  true  earth  radius 

is  a  repr  esentative  value  for  the  surface  index  of  refraction 
and  may  be  consider  ed  to  be  unity  for  all  pr  actical  purposes 
is  the  effective  gradient  of  radio  refractivity  with  height 


,  ,  N 

where  n  =  1  +  and 


is  the  local  elevation  angle  of  the  ray 


In  general  the  vertical  index  gradient  is  not  independent  of  height, 
in  which  case  Ae  would  have  little  analytical  usefulness  However,  for  the 
moment,  let  an  equivalent  effective  earth  radius  Ae  be  defined  to  satisfy  a 
particular  propagation  condition  observed  from  the  interferometer  results 


Referring  to  Figure  1  and  Equation  (3),  at  a  particular  time  when  a 
fade  occurs  there  is  only  one  value  for  the  earth  radius  which  satisfies  the 
measurement  of  IDn,  h1l  and  hg'  to  produce  the  n**1  fade  Therefore,  at  each 
fade  condition  the  geometrical  measui  ements  will  generate  a  unique  value 
for  Ae  which  will  generally  be  different,  from  fade-to-£ade  The  only  time 
the  values  for  Ae  could  be  constant  is  A  one  value  for  Ae  can  be  found  which 
satisfies  the  geometry  for  all  ranges  and  initial  elevation  angles  This 
situation  occurs  only  if  the  vertical  gradient  dN/dh  is  everywhere  a  constant, 
a  condition  not  generally  representative  in  reality,  and  where  cos  0  is  taken 
as  unity  for  propagation  at  small  elevation  angles  (Equation  (4)) 


6 


3. 


A  COMPARISON  OF  INTERFEROMETER  AND  RAY  TRACING 
CALCULATIONS  OF  THE  EQUIVALENT  EFFECTIVE  EARTH'S  RADII 


A  total  of  eighty  pattern  measurements  was  carried  out  at  the  Royal 
Canadian  Air  Force,  Primrose  Lake  Evaluation  Range,  Cold  Lake,  Alberta, 
Canada.  2  The  receiving  site,  using  a  wide-angle  microwave  antenna,  was 
located  on  a  hill  overlooking  a  large  fresh  water  lake  The  radio  source  was 
carried  in  an  aircraft  at  constant  altitude  (as  well  as  could  be  flown)  with 
theodolite  cameras  and  radar  altimeter  recordings  used  for  the  final  accurate 
height  determinations.  The  electrical  length  of  the  direct  path  between  the 
receiver  and  the  aircraft  was  measured  by  a  short  pulse  radar  and  time 
interval  electronics  system  located  at  the  receiving  site. 

Independent  radio  refractivity  measurements  were  made  using  an 
instrumented  aircraft  together  with  periodic  radiosonde  soundings  during  the 
time  the  radio  experiment  was  being  conducted.  The  refractivity  data  were 
then  used  in  a  ray  tracing  analysis  to  determine  the  angles  of  arrival  and 
electrical  path  lengths  corresponding  to  the  position  of  the  signal  source  air¬ 
craft  for  each  radio  fade  condition  The  ray  tracing  results  were  then  used 
to  further  calculate  "equivalent  effective  earth  radii"  which  would  satisfy  the 
initial  elevation  angles  and  ranges,  corresponding  to  the  geometry  at  each 
fade 


The  radio  propagation  and  independent  ray  tracing  results  for  Ae  were 
then  compared  for  all  the  tests  undertaken. 

Figure  4  shows  the  average  differences  between  the  results,  together 
with  the  standard  deviation,  plotted  against  the  transmitter  range  R.  The 
estimated  error  in  the  interferometer  results  is  shown  by  the  dashed  curve 
(Appendix  II)  The  errors  increase  at  short  ranges  due  to  the  fact  that  the 
signal  source  height  and  range  errors  become  more  significant. 

One  objective  in  the  experiment  was  to  relate  the  radiometric  and  ray 
tracing  calculations  of  the  equivalent  effective  earth  radii,  Ae,  to  within 
200  nm.  On  the  average  these  results  were  nearly  achieved  and  to  do  this 
required  a  high  degree  of  accuracy  such  as  signal  source  height,  h2,  within 
±25  feet,  receiver  height,  hlf  within  ±0  05  feet,  signal  source  range,  R, 
within  ±100  feet  and  a  wavelength  determination,  \,  within  ±0.001  feet. 
During  the  experiment  both  vertical  and  horizontal  polarizations  were 
available  for  use  at  L-  and  S-Bands 
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FIGURE  4  EFFECTIVE  EARTH  RADIUS  DIFFERENCE  (nm)  BETWEEN 

RAY- TRACING  AND  INTERFEROMETER 
RESULTS  VERSUS  RANGE 


4. 


A  COMPARISON  OF  EFFECTIVE  PROPAGATION  CONDITIONS  WITH 
RADIOMETRIC  MEASUREMENTS 


Any  one  set  of  earth  radii  values  generate  a  corresponding  set  of  re¬ 
fraction  gradients,  dN/dh,  as  shown  from  Equation  (4)  Larger  values  of 
earth  radii  represent  correspondingly  larger  values  for  the  magnitudes  of 
the  gradients. 

Figure  5  is  a  plot  of  the  received  radio  signal  levels  for  a  particular 
flight  Figure  6  shows  the  radiometric  calculations  for  Ae,  obtained  for  each 
fade  condition  plotted  against  the  radio  source  range  from  the  receiving  site. 
Two  sets  of  data  are  shown  corresponding  to  inbound  (Figure  5)  and  outbound 
flights  Within  thirty  miles  range  the  accuracy  decreases  due  to  the  effect  of 
range  and  height  errors,  as  previously  discussed  Beyond  thirty  miles  the 
magnitudes  of  the  earth  radii  values,  Ae,  are  seen  to  be  essentially  constant, 
which  suggests  that  a  constant  gradient  refractivity  profile  represents  the 
propagation  conditions.  Figure  7  shows  the  radio  refractivity  profile  measured 
with  a  radiosonde  from  a  launch  point  approximately  twenty-five  miles  from 
the  receiving  site  and  in  the  direction  of  the  aircraft  flight  path.  As  indicated 
the  refractivity  gradient  is  essentially  constant. 

Figure  8  shows  a  different  situation,  where  the  radiometric,  Ae, 
values  are  seen  to  increase  with  range,  which,  therefore,  suggests  that  the 
gradient  is  increasing  in  magnitude.  Figure  9  shows  the  measured  refrac¬ 
tivity  profile  and  it  is  apparent  that  the  gradient  is  increasing  at  the  lower 
levels.  As  the  signal  source  range  increases  the  ray  passes  through  more  of 
the  lower  atmosphere,  where  the  larger  gradient  would  become  more  effective 
Therefore,  the  radiometric  data  and  the  radio  refi  activity  profile  characters 
istics  show  a  consistent  behavior 

In  general,  the  radiometric  and  propagation  data  show  surprisingly 
good  correspondence  in  terms  of  the  relative  variations  of  earth  radii  and 
the  effective  refractivity  gradients  with  height.  Aircraft  meteorological 
measurements  were  similar  to  the  radiosonde  data  with  the  exception  that 
greater  fine- structure  could  be  resolved  This  latter  information  is  not 
particularly  important  to  the  average  conditions  under  discussion  herein. 

The  remaining  problem  is  how  to  exploit  the  radiometric  measurement 
technique  to  define  the  effective  propagation  conditions  Several  possible 
configurations  are  discussed  in  the  following  sections. 
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FIGURE  5  SIGNAL  AMPLITUDE  VERSUS  RANGE 
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FIGURE  8  EFFECTIVE  EARTH  RADII  VALUES  VERSUS  RANGE 
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5. 


A  DISCUSSION  OF  METHODS  TO  USE  RADIOMETRIC  MEASURE¬ 
MENTS  TO  DETERMINE  PROPAGATION  CONDITIONS 


Any  determination  of  propagation  conditions  is,  in  the  first  place,  re¬ 
stricted  to  the  maximum  height  of  the  signal  source  aircraft  However, 
since  most  of  the  non-standard  variations  in  the  refractivity  profile  occur 
below  15,  000  feet  the  experiment  can  be  limited  in  height;  with  radiosonde 
measurements  used  above  this  level  to  define  the  average,  spherically 
stratified  conditions 

In  order  to  reduce  the  effect  of  range  measurement  errors  and  the 
cos  0  factor  (Equation  (4)),  the  source  should  be  located  at  long  ranges  but 
within  the  radar  horizon  The  larger  the  wavelength  the  greater  will  be  the 
separation  between  fades  However,  the  number  of  fades  will  be  correspond 
mgly  reduced  which  limits  the  precise  definition  of  the  propagation  conditions 

As  the  aircraft  range  changes,  and  with  the  height  maintained  con¬ 
stant  the  degree  of  ray  path  curvature  depends  upon  the  initial  elevation 
angle  which,  in  this  case,  is  the  angle  of  arrival  In  order  to  obtain  a 
meaningful  measurement  of  the  propagation  conditions  it  would  be  convenient 
to  generate  the  values  of  effective  earth  radii  as  a  function  of  radio  source 
height  rather  than  in  terms  of  range. 

These  considerations  suggest  an  experiment  where  the  radio  source 
is  constrained  at  some  desirable  range  and  rises  from  a  point  below  the 
radio  horizon  to  some  maximum  altitude 


As  the  source  rises  above  the  radio  horizon  the  first  order  fade 
(n  =  1)  can  be  determined  and  all  subsequent  order  fades  then  accountable. 


From  Equation  (3)  fades  occur  whenever 
,  nRX 


ai*d 


Zh^ 

n  =  1 ,  2,  3,  .  ,  .  ,  etc . 


Replacing  h2‘  and  hx!  by  their  c or  r  esponding  values  of  hx  and  h2 
(Appendix  I),  gives 


L  £ 


(R  -  nil 

2Ae 


nX  R 


2  ^  -  It 


(5) 


(6) 
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and 


(R  -  r ! )  (hi  - 


(?) 


ri  ~ 


h2 


_ 2Ae 

(R  -  EZ 

2Ae 


which  for  a  particular  value  of  n,  X;  R,  hlf  and  h2,  generates  an  earth  radius 
Ae 


In  the  Cold  Lake  analyses,  equations  of  this  kind  were  solved  using  a 
computer  due  to  their  complexity  One  may  obtain  an  estimate  of  the  varia¬ 
tion  of  received  signal  with  height  using  Equation  (5)  and  making  the  tacit 
assumption  that  h^1  and  h2’  are  sufficiently  close  to  hx  and  h2  to  let 


h2  ~ 


nRX 

Zhx 


(8) 


For  example,  with 


R  =  100  km 
X  =1  meter 
hx  -  50  meters 

the  fades  will  occur  for  source  heights  of  1000  meter  intervals 

With  the  source  at  1000  meters,  Equation  (8)  shows  that  the  aircraft 
range  if  variable,  would  produce  a  first  order  fade  at  ZOO  km,  second  fade 
at  100  km  and  a  third  fade  at  66  km  At  gi  eater  source  heights  the  range 
variations  required  to  produce  range  iependent  fades  will  require  even 
greater  range  changes.  Therefore  it  would  be  relatively  easy  to  produce 
a  height-dependent  fade  pattern  without  placing  severe  constraints  on  the 
permissible  aircraft  range  variations  The  height  intervals  which  produce 
fade  data  can  be  decreased  of  course,  by  raising  the  receiver  height  or  by 
decreasing  the  wavelength  However  these  changes  will  affect  the  per¬ 
missible  variation  of  range  R  and  may  not  be  acceptable  in  practice 

The  remaining  problem  is  to  use  the  recorded  and  calculated  data  to 
determine  an  effective  refractivity  profile.  Referring  to  Figure  10,  the 
following  procedure  could  be  attempted  The  source  rises  above  the  radio 
horizon  producing  the  first  fade  (n  =  1)  at  which  time  the  range  R,  and  air¬ 
craft  height,  (h2)0l  is  recorded  The  solution  of  the  geometry,  (Equations  (6) 
and  (7)),  gives  a  particular  value  for  the  effective  earth  radius,  Ae  From 
Equation  (4)  with  n0  and  cos  0  taken  as  unity  an  effective  gradient  (dN/dh)01 
is  determined  This  gradient  will  then  represent  the  average  propagation 
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FIGURE  10  PIECEWISE  CONSTRUCTION  OF  THE  RADIO 
REFRACTIVITY  PROFILE 
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condition  over  this  first  height  interval  (h2)01  and  over  the  horizontal  extent 
of  the  range  R  When  the  second  fade  occurs  (n  =  2),  in  the  same  way  the 
average  effective  gradient  (dN/dh)02is  found  which  extends  over  the  height 
interval  (h2)02  . 


It  is  readily  seen  that 


dN 
(  dh 


OP 


dN 

(h2)op^  (— )  (ha)< 

dh  01 


dN  v 
11  +  (^h  1 


(hs) 


IS 


IS 


(9) 


from  which  the  average  effective  gradient  over  the  interval  (h2)l2canbe 
determined  The  procedure  is  continued  with  each  fade  in  height  until  the 
overall  vertical  gradient  profile  is  determined  in  height  sections  Also,  the 
relative  values  for  the  radio  refractivity  can  then  be  determined  from  a 
knowledge  of  the  individual  gradients  and  the  magnitudes  of  the  respective 
height  intervals.  In  theory,  a  vertical  refractivity  profile  can  then  be  con¬ 
structed  and  made  to  fit  the  appropriate  refractivity  readings  obtained  at 
higher  elevations  with  the  radiosonde  data  Since  the  surface  refractive 
index,  n0,  is  probably  quite  variable  it  would  be  difficult  to  select  a  suitable 
n0  from  which  to  generate  the  absolute  profile  For  this  reason  it  is  more 
meaningful  to  tie  the  profile  to  measurements  at  the  higher  levels.  These 
high-level  measurements  can  be  obtained  from  a  radiosonde  launched  any¬ 
where  within  the  general  area. 
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6 . 


THE  SENSITIVITY  OF  THE  RADIOMETRIC  MEASUREMENTS  TO 
CHANGES  IN  THE  REFRACTIVIT Y  GRADIENTS 


The  results  obtained  in  the  Cold  Lake  experiment  demonstrated  that 
there  was  a  general  tendency  for  the  effective  earth  radii  values,  Ae,  to  vary 
in  accordance  with  changes  in  the  propagation  conditions  It  is  useful  to 
analyze  the  expected  variations  as  a  function  of  gradient  changes  in  order  to 
determine  whether  or  not  sensible  data  may  be  obtained  within  limits  imposed 
by  experimental  errors 

Assuming  that:  the  same  accuracies  can  be  obtained  as  in  the  Cold 
Lake  experiment  (Section  3),  the  error  in  the  effective  earth  radii  calcula¬ 
tions  will  be  set  at  +200  nautical  miles 


From  Equation  (4),  with  cos  0  set  equal  to  unity,  then 


6(Ae)  Ae  ,  dN 

~aT~  ~  '  To3"  6,'dh) 


(10) 


where  N  is  the  radio  refractivity 

Ae  is  the  effective  earth  radius 


Considering  standard  conditions,  then  setting  Ae  to  be  4600  nm,  (i.  e.  , 
8500  km)  then 

6  (dN/dh)  =  +5N/km  (11) 

For  experimental  conditions  outlined  in  Section  5,  this  result  indicates 
that  the  gradient  calculation  over  the  first  kilometer  height  interval  could  not 
be  determined  to  better  than  +5N/km.  Comparing  with  an  exponential  refer¬ 
ence  atmosphere,  3  where  the  initial  gradient  is  about  40  N/km  near  sea  level, 
it  is  apparent  that  the  measurement  would  be  in  error  by  at  least  ±  12  percent. 

At  the  second  height  interval  the  average  gradient  error  would  again 
be  about  12  percent.  Referring  to  Equation  (9)  the  gradient  in  the  interval 

would  be  affected  by  the  combined  err  ors,  and  so  on>  with  the  cumulative 
errors  increasing  with  height.  Since  there  is  no  bias  toward  the  sign  of  the 
errors  the  final  accumulated  error  over  five  fade  conditions  would  not  be 
expected  to  be  60  percent. 
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7. 


SOME  COMMENTS  ON  THE  ERROR  MAGNITUDES  AND  METHODS 
TO  REDUCE  THEM 


From  the  error  analysis  equations  presented  in  Appendix  II,  the  Cold 
Lake  data  clearly  show  that  the  inability  to  determine  the  signal  source  height 
contributed  to  most  of  the  error. 

The  development  of  radar  altimetry  provides  a  present  day  capability 
to  measure  the  height  above  a  smooth  water  surface  to  within  a  few  feet.  If 
the  radar  is  placed  on  a  stabilized  mount  to  ensure  that  the  range  measure¬ 
ment  is  in  the  vertical  an  experiment  could  be  conducted  without  imposing 
severe  flight  requirements  on  the  aircraft  The  effect  of  sea- surface  rough¬ 
ness  should  average  out  since  the  radar  antenna  pattern  will  encompass  a 
significantly  large  area  on  the  sea  surface 

With  a  receiver  height  hx  of  50  meters,  for  example,  the  reflection 
area  or  first  Fresnel  Zone  will  be  very  large  and  located  several  miles  from 
the  shore  line  This  has  two  advantages  in  that  surface  roughness  is  aver¬ 
aged  out  and  the  effect  of  local  upwelling  of  the  water,  near  the  shore  line, 
is  reduced 

The  difficult  problem  is  to  determine  the  height,  hv  of  the  receiver, 
relative  to  the  reflection  area  One  method  to  derive  this  effective  height 
would  be  to  transmit  two  signal  frequencies  from  the  aircraft  and  to  compare 
the  heights  at  which  the  first  order  fades  occur  for  the  individual  frequencies 

Then,  using  the  Equations  (6)  and  (7)  the  first  order  fades  occur  where 

^  (Rx»  hi »  Ae)  (1^) 

^  (R  s?  ^  i,  Ae)  (13) 

where  the  difference  between  and  X2  is  only  large  enough  to  separate  the 
individual  signals  at  the  receiver.  If  the  fades  occur  within  a  small  variation 
of  range,  R,  and  height,  h2,  then  Ae  can  be  considered  to  be  the  same  for 
both  conditions  The  effective  receiver  height,  hlf  can  then  be  found 

There  may  be  some  advantage  to  continue  the  use  of  two  frequencies 
since  several  values  for  hx  can  then  be  calculated  and  a  more  meaningful 
average  value  determined  Furthermore,  the  use  of  two  frequencies  pro¬ 
vides  a  second  measure  of  the  refractivity  profile  at  different  heights  than 
the  values  obtained  with  the  first  frequency  In  the  event  of  antenna  shading 
by  the  aircraft  the  second  set  of  data  provides  additional  reliability  in  the 
experiment . 
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However,  from  the  error  analysis  Equation  (3)  (Appendix  II)  assuming 
two  wavelengths  are  used,  and  with  R  and  Ae  essentially  the  same  for  both 
fade  conditions  at  Xa  and  Xs  then 


u  R  (^g ~ 

2  ((h8)8  -  (haJj) 


(14) 


For  a  range  of  100  km,  a  wavelength  separation  of  0  1  meters  at  a 
mean  wavelength  of  1  meter,  then  estimating  from  Equation  (8)  the  fades 
occur  at  1000  and  1100  meters  Assuming  a  range  measurement  error  of 
4-30  meters,  a  wavelength  measurement  error  of  4-0  0003  meters  and  a  height 
measurement  error  of  ±  2  meters,  then  the  errors  in  hi  from  Equation  (14)  are 


a  0  015  m  due  to  the  range  uncertainty, 

b  0,3  m  due  to  the  wavelength  uncertainty, 

c  10m  due  to  the  aircraft  height  uncertainty 


or  an  accumulated  error  of  about  1  3  meters  largely  due  to  the  aircraft 
height  error 


From  Equation  (11)  on  Appendix  II,  this  leads  to  an  estimate  of  the 
effective  earth  radius  error  due  to  ha  where 


6  Ae 
(  Ae 


2haAe  6hi 

=  11  ■  '1T> 

=  (0.7)  <6h,/ht) 

=  (0  7)  (1/50)  X  100 


=  1  8%  er  r  or 


(15) 


This  is  an  acceptable  error  and  it  is  seen  that  the  effect  of  the  receiv¬ 
ing  height  uncertainty  is  reduced  under  long  range  propagation  conditions  and 
by  placing  the  receiver  at  a  significant  height  above  the  reflecting  surface. 

From  the  error  analysis,  Equations  (10),  (11),  ( 1 Z ) P  (13),  (Appendix 
II),  the  error  contributions  to  the  determination  of  the  effective  earth  radii 
are  then 


a 

0  11% 

due 

to 

a  30  meter  range  error. 

b. 

1  8% 

due 

to 

a  1  3  meter  receiving  antenna  height  error 

c  , 

0  34% 

due 

to 

a  2  meter  aircraft  height  error 

d. 

0.  051% 

due 

to 

a  0  0003  meter  wavelength  error. 

That  is,  a  total  error  of  2.  3% 
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From  Equation  (10)  the  gradient  error  is  then  0  27%  which  for  the  ex¬ 
ponential  reference  atmosphere  represents  a  gradient  error  of  about  ±0.  11 
N  units  per  km<  This  accuracy  is  certainly  adequate 

It  is  significant  to  note  that  the  bulk  of  this  error  is  contributed  by  the 
indeterminacy  in  the  receiver  height  relative  to  the  effective  reflecting  sur¬ 
face.  Over  a  smooth  water  surface,  such  as  an  inland  lake,  the  receiver 
height  can  be  determined  very  accurately  through  the  use  of  standard  survey¬ 
ing  methods  However,  over  the  ocean,  the  height  would  probably  have  to  be 
indirectly  determined  such  as  with  the  two  frequency  method  in  which  case 
the  receiver  height  error  is  expected  to  remain  the  more  difficult  parameter 
to  calculate. 


From  the  error  analysis  Equation  (11)  (Appendix  II)  it  is  seen  that 


6Ae  2h?Ae  .  6h , 

'ir’h,5  11 


(16) 


At  the  first  order  fade,  (n  =  1),  h2  equals  1  km,  range  R]  equals  100 
km,  then 


2hg Ae  _ 

— zr? >  1 

R2 


Therefore,  as  successive  data  are  obtamed  at  greater  heights,  the 
receiver  height  error  6h1?  will  become  larger  Therefore,  it  is  of  particular 
importance  that  hx  be  determined  by  any  possible  method  which  will  improve 
its  accuracy 

Furthermore,  the  gradient  decreases  with  height  (based  upon  an  ex¬ 
ponential  model)  in  which  case  the  above  error,  Sh^  becomes  even  more 
significant  Around  5  km,  in  height  a  typical  refractivity  variation  is  20  N 
units  per  kilometer. 

From  the  error  analysis  and  using  Equation  (15)  the  accumulated  error 
after  five  fades  is  readily  shown  to  be  about  1  4%  in  the  effective  gradient 
determination  for  the  5  kilometer  height  section  (h2)06.  Similar  to  the  method 
illustrated  by  Equation  (9)  the  error  in  the  gradient  determination  over  the 
height  interval  (hs)^  is  then  given  by  the  accumulated  errors  in  (dN/dh)  over 
each  section  and  weighted  by  the  factor,  h^,  as  shown  by  Equation  (16).  The 
accumulated  error  is  readily  shown  to  be  about  4%,  which  is  acceptable . 
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SUMMARY  AND  CONCLUSIONS 


The  vertical  reflection  interferometer  experiment  conducted  at  Cold 
Lake  demonstrated  that  radiometric  measurements  could  provide  values  of 
the  effective  earth  radii  which  compare  favorably  with  independent  ray  trac¬ 
ing  analyses. 

These  results  lead  to  a  procedure  which  could  generate  the  effective 
vertical  radio  refractivity  profile  by  a  modified  experimental  procedure  In 
this  case  a  radio  source  is  carried  in  height  at  a  long  range  from  the  receiver 
such  as  to  generate  a  multipath  fading  pattern  which  is  essentially  only  height 
dependent.  At  the  time  of  each  ordered  fade,  the  solution  of  the  geometry 
generates  a  particular  earth  radius  value,  Ae  which  can,  in  turn,  be  used 
to  determine  an  effective  refractivity  gradient  over  the  corresponding  height 
interval  Successive  calculations  at  each  fade  condition  lead  to  a  piecewise 
construction  of  the  overall  refractivity  profile, 

In  order  to  improve  the  experimental  accuracy  the  signal  source 
height  relative  to  the  reflecting  surface  should  be  determined  with  the  most 
accurate  radar  altimeter  equipment  A  second  difficulty  which  would  be 
experienced  in  practice  is  to  determine  the  effective  height  of  the  receiving 
antenna  relative  to  the  reflecting  surface.  It  is  shown  that  it  can  be  ac¬ 
complished  by  observing  the  relative  fading  patterns  generated  by  using  two 
simultaneous  frequencies  In  fact,  under  long  range  conditions  and  with  a 
receiver  placed  well  above  the  reflecting  surface,  the  errors  in  this  cal¬ 
culation  would  be  acceptable 

The  effects  of  anomalous  propagation  conditions  were  not  considered 
since  it  is  apparent  that  the  experiment  would  break  down  should  either  the 
direct  or  reflected  signal  be  removed 
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APPENDIX  I 

THE  "EFFECTIVE  EARTH  RADIUS"  CONCEPT 
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THE  "EFFECTIVE  EARTH  RADIUS"  CONCEPT 


Referring  to  Figure  11,  a  receiver,  ’A',  at  height  hx] ,  will  receive 
signals  along  two  paths  from  a  transmitter,  ’F1,  at  height  h2’  The  first  is 
the  direct  path  F  to  A,  and  assuming  that  the  surface  is  highly  reflective, 
the  second  path  is  FGA  As  the  transmitter  moves  in  range,  the  signals 
combine  at  the  receiver  to  produce  minima  and  maxima  amplitude  variations 
For  a  water  surface  and  with  small  grazing  angle,  Bf  the  phase  of  the  re¬ 
flected  signal  is  changed  by  1  80  degr  ees,  i .  e  ,  radians.  The  fir  st  fade  in¬ 
side  the  radio  horizon  is  given  the  order  n  =  1  and  it  occurs  at  a  correspond¬ 
ing  transmitter  range  D  v  As  the  transmitter  moves  towards  the  receiver, 
each  fade  condition  (minimum  received  signal)  occurs  for  ranges  given  by 


D 


n 


zy  h3' 

nX 


U-1) 


where  n  =  1,  2,  3,  .  .  ,  ,  etc  „ 

X  =  the  radio  signal  wavelength 


In  the  real  situation,  of  course,  the  reflector  becomes  the  surface  of 
a  spherical  earth  As  shown  in  Figure  12,  the  true  heights  of  the  receiver 
and  the  transmitter  over  the  spherical  surface  are  AA1  and  BB',  respectively 
If  a  plane  is  passed  through  the  reflection  point  C,  Equation  41— J.)  can  be  used 
to  determine  the  fade  dependence  upon  range  if  h and  h?  are  expressed  in 
terms  of  hl  and  respectively. 


From  triangle  BCO, 

Q  ( A  +  h  g ) 2  ••  A s  -  r?2 

sm  p  =  - - — - - 


also  from  triangle  BCB" 
h?!  -  r2  sin  fi 


Substituting  for  sin  into  Equation  (1-3)  and  simplifying  gives 


h 


\ 

3 


(1-2) 


(1-3) 


(1-4) 


Over  the  range  of  experimental  interest  r?will  be  very  much  greater  than  h5 
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FIGURE  II  FLAT  EARTH  REFLECTION  GEOMETRY 
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FIGUF  12  GEOMETRY  OF  FLAT  EARTH  APPROXIMATION 
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Therefore 


hs'  ~  hs  - 


2A 


(1-5) 


similarily  it  can  be  shown  that 

S 


hi'  »  «i 


hi  -  ^ 


2  A 


d-6) 


Substituting  these  equations  into  Equation  (1-1)  gives 


Dn  a  2  (h>  '  *£)  ^  ^ 


d-7) 


Since  the  angles  of  incidence  and  reflection  are  equal,  we  have 


hi  - 


r  l 


ZA 


hp  - 


rs 


d-8) 


ZA 


However,  at  any  particular  fade,  of  order  n,  we  have 
(r x  +  r a)  -  Dn  =  nX 
or  r?  =  Dn  -  rx  +  nX 


(1-9) 

(1-10) 


Over  the  range  of  experimental  interest,  nX  will  be  less  than  about 
ZO  feet  (n  <  50  for  X~10  cm  and  n  <  Z0  for  X— 30  cm)  and  rl  will  be  greater 
than  10  miles  Therefore  Equation  (1-10)  is  approximately 


rs~  Dn  -  rl 

Substituting  into  Equation  (1-7)  and  (1-8)  gives 


(I- 1  1) 


D 


n 


4 


nX 


and 


D 


= 


n 


ZA 


hs  - 


jDn  Z-llL 

2A 


(1-12) 


(1-13) 
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Finally  from  triangle  AAnC  and  Equation  (1-6) 


sin 


d-14) 


Under  propagation  conditions  in  the  real  atmosphere,  of  course,  the 
range  Dn  to  the  transmitter  is  a  curved  path  rather  than  a  straight  line  How¬ 
ever,  if  a  straight  path  range  Dn  is  assumed,  and  ihe  terminal  heights  hl  and 
h3  retained  in  magnitude,  the  geometry  as  shown  in  Figure  12  will  simply 
give  a  value  of  the  earth's  radius  Ae  which  is  larger  than  the  true  earth 
radius  A  The  above  equations  are  therefore  quite  general  and,  by  substitut¬ 
ing  Ae  for  A,  they  can  be  used  to  generate  values  of  the  effective  earth's 
radii  x  A  7094  computer  program  was  used  to  solve  the  above  equations  for 
the  particular,  known  input,  parameter  s  X,  n,  Dn,  hx,  and  h2  involved  in  this 
experiment  From  the  resulting  tables  the  appropriate  values  for  Ae  can  be 
selected,  corresponding  to  the  exper  mental  conditions  at  any  particular  fade 
of  order  n„ 
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INTERFEROMETER  ERROR  ANALYSIS 
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INTERFEROMETER  ERROR  ANALYSIS 


The  error  introduced  in  the  Ae  calculations  can  be  determined  based 
on  the  experimental  accuracies  oi  Dn,  h1#  and  h?.  To  simplify  this  analysis* 
an  approximate  relationship  between  these  parametei  s  can  be  detei  mined 
from  the  geometry  of  Figure  13.  This  geometry  was  not  used  for  Ae  calcula¬ 
tions  (Reference  Appendix  I)  because  it  is  slightly  less  accurate  and  does  not 
give  all  the  parameters  required  in  the  overall  analysis  such  as  r1  and  rP „ 
Under  the  range  of  experimental  conditions  of  interest  it  may  be  assumed 
that  the  reflection  point  is  very  near  the  receiving  terminal  and  the  earth  is 
essentially  flat  within  this  confined  area  Therefore  the  flat  earth  Equation 
(I- 1 )  can  be  applied  From  Figure  13,  it  can  be  shown  that 


h Pf  -  h? 


2 A  4"  hp 
2  (A  +  hA) 


D 


n 


U  2 


2  (A  +  hx)  2  (A  +  hx) 


(II-l) 


and  for  hx  <<:A  and  hp  <<A  Equation  (II - 1 )  becomes 


hp1 


D 


n 


2A 


(II-2) 


Substituting  this  equation  into  ^1-1),  which  becomes 


D 


n  ~ 


nX 


(since  h1!  =  h1),  we  find  that 
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Again,  this  equation  is  general  and  A  can  be  replaced  by  Ae  to  obtain  the 
dependence  of  Ae  on  the  parameters  measured  and  applied  in  the  interfero¬ 
meter  analysis.  Then, 


log  Ae  =  2  log  Dn  +  log  hL  -  log  (2hxhp  -  n>  Dn)  (II-4) 

The  effect  of  errors  in  each  of  the  separate  parameters  Dn,  h1#  hp,  and  X  can 
be  found  by  differentiating  Equation  (II-4)  to  give 
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From  Equation  (11-3) 
2hih2  -  nXDn  ~ 
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(II- 9) 


Substituting  this  equation  into  the  above  gives 
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Table  II- 1  shows  the  individual  errors  in  6Ae/Ae  for  various  test 
conditions  and  the  total  error  equivalent  to  the  sum  of  the  individual  errors 
Finally,  multiplying  this  total  error  by  average  earth  radii  conditions  gives 
the  estimated  maximum  errors  in  the  Ae  calculations  A  graph  of  this  error 
relationship  to  transmitter  range  is  shown  in  Figure  4 
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